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ENERGi  EQUATION  APPROXIMATIONS  IN  STC2D  MECHANICS 
By  Arthur  W.  Goldstein 
Levis  Research  Center 

National  Aeronautics  and  Space  Administration 
Cleveland,  Ohio 

There  is  some  confusion  in  the  literature  of  fluid  mechanics-*-*^  in  re¬ 
gard  to  the  correct  form  of  the  energy  equation  for  the  study  of  the  flow 
of  nearly  incompressible  fluids .  Several  forms  of  the  energy  equation  and 
their  use  are  therefore  discussed  in  this  note. 

The  general  energy  equation  for  a  fluid,  when  the  chemical  composition 
is  fixed,  may  be  expressed  in  terms  of  the  entropy  S,  the  enthalpy  H,  or 
the  internal  energy  C  as 

pTS  =  pH  -  p  =  pH  +  ppv  =  J  (l) 


i 


where 

p  =  l/v,  density 
T  =  temperature 
p  =  pressure 

<7  -  net  rate  of  heat  addition  per  unit  volume,  including  conduction  (-V’kVT), 
dissipation  (t),  heat  sources  (pQ,) 

The  dot  indicates  the  substantial  time  derivative,  e.g., 

U  a  dU/dt  =  dU/cSt  +  V  *  VC 

When  the  energy  equation  is  used  in  this  form,  no  question  arises,  but  when 
H  or  U  is  expressed  in  terms  of  T,  an  error  can  appear  by  neglecting  the 
terms  (dH/dp)„,  or  (dr/dv).,,. 
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The  commonly  used  coefficients  needed  for  the  discussion  are: 


Cp  2  (c*H,/c)ll)p  =*  T(dS/dT)p,  specific  heat,  constant  pressure 
Cv  =  (du/dT)  a  T(dS/dT;v-,  specific  heat,  constant  volume 
a  =  (dv/dT)  /v,  coefficient,  of  t-hennal  expansion 

Jr 

=  '•(dv/dp)fn/vj  coefficient  of  compressibility 
j3  a  (dp/dT)y/p,  coefficient  of  tension 


The  existence  of  an  equation  of  state  between  p,  v,  T  implies 


a 

»  PP-jP  (Eq. 

(1.05))3 

(2) 

The  reci 

.procity  relation  is 

(si),. 

■  (§f)Y  ■ » 

(Eq.  (4.23))3 

(3) 

or 

(§,- 

-(§$  =  -av 
\  /  p 

(Eq.  (4.46))3 

(4) 

Also 

Gp  - 

Gy  =  aSpvT 

(Eq.  ( 4.27 ) )3 

(5) 

If  H  = 

H(T,p)  and  “j  «  *-v~. 

,r),  t-hen 

H  =■  CpT  +  (l  - 

•  aT)vp 

(6) 

i’  a  C-yT  -  (l  - 

-  0T)pv 

(7) 

and  with  the  above  definitions,  the  energy  equation  may  be  stated  in  the  two 


forms 


v  «  pCv*  +  PTppv 
«T  *»  pCpT  -  aTp 


For  a  perfect  gas,  a~  »  J3T  «  1,  yielding  the  familiar  forms , 


H  -  CpT,  U  -  CVT,  »T  -  pCpT  -  p 


PC- 


m 


ppv 


(8a) 

(8b) 
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(9) 
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With  liquids  it  would  seem  reasonable  to  neglect  the  terms  involving  v 
in  Eqs.  (7)  and  (8a).  This,  however,  is  not  generally  permissible  because 
3  can  become  very  lar&e.  For  liquids  of  very  small  compressibility  ({^ 
small),  the  value  of  p(  =  a/p^p)  may  be  large  with  a  consequence  that  the 
term  pTppv  in  (8a)  and  PTpv  in  (7)  might  not  be  negligibly  small.  The 
kind  of  magnitudes  involved  may  be  Judged  from  the  data  of  the  table: 


Liquid 

i - 

P; 

a, 

°C 

atm 

1 

°c 

*2° 

220 

50 

0.001546 

Hg 

100 

1.0 

0.000180 

gm  I  gm-°C  I  °C 


The  table  would  suggest  that  for  many  liquids  (3T  »  1.  Eq.  (7)  may 
then  be  expressed  as  U  »  CyT  +  PTpv  =  pTS  =  J  to  the  effect  that  all  heat 
addition  increases  internal  energy  and  none  is  available  for  external  work. 


The  energy  equations ( 8a, b)  do  not  permit  such  general  simplification  inde¬ 
pendently  of  the  particular  problem*  For  example,  in  a  problem  where  p/p 
is  of  the  same  order  as  T/T,  the  ratio  of  the  term  aTp  to  pCpT  in 
Eq.  (Gb)  is  of  the  order  of  ap/pCp,  which  for  water  and  mercury  at  condi¬ 
tions  cited  in  the  table  has  values  of  2X10"3  and  10"^,  respectively.  Por 
these  fluids  under  such  circumstances,  Eq  (8b)  becomes 


/rt 


j  «  pCjf 


(10) 


On  the  other  hand,  for.  isentropic  flows  (such  as  an  accoustic  wave),  J  is 

A 

small  relative  to  both  pCpT  and  aTp,  and  both  these  terms  must  be  retained* 
In  many  flow  problems,  the  use  of  Eq.  (8b)  is  preferable  to  (8a)  because 


the  ranges  of  p  and  T  are  indicated  by  the  boundary  conditions,  making 
it  possible  to  decide  in  advance  when  either  term  may  be  omitted* 


It  is  often  assumed  for  liquids  that  Cp  Cv  (e.g.,  Ref.  3,  p.  4l) 
because  of  the  factor  av  in  Eq.  (5).  The  effect  of  the  additional  factor 
PpT  may  render  this  assumption  inadequate  in  many  cases  as  shown  in  the 
table.  In  cases  where  (8b)  may  be  simplified  to  (10),  the  corresponding 
simplification  of  (8a),  namely  J  =•  C^p$,  is  not  admissible  unless  the 
liquid  satisfies  the  additional  condition  Cp  »  cv>  which  cannot  be  assumed 
for  all  condensed  fluids.  The  particular  fluid  and  problem  together  will 
determine  whether  (8a)  or  (8b)  yields  the  simplest  equation. 

The  approximation  used  by  Millsaps  and  Pohlhausen1  is  J.  =  pCvT  for 
water  at  30°  C  and  unstated  pressure,  which  will  be  assumed  to  be  1  atm. 

For  these  conditions,  a  =  0.0003/°C,  (3^  =  48xi0“6/atm,  from  which 
(Cp  -  Cy)/Cp  =  0.014.  Since  the  effect  of  p  is  small,  the  correct 
approximation  is  J  =  pCpT.  The  discrepancy  in  the  energy  equation  is 
minor  because  of  the  small  value  of  (Cp  -  Cv)/Cp.  The  same  equation  is 
also  used  for  argon  where  the  error  is  not  small,  since  (Cp  -  Cy)/Cp  =  0.4. 
The  proper  approximation  with  small  p  is  J  *»  pCp$  as  was  noted  by 
Ostrach  and  Thornton.5 

Hidalgo2  uses  the  approximation  J  =  pCpT  -  £  for  liquid  glass.  In 
principle,  the  equation  is  not  correct,  since,  if  the  p  term  is  to  be 
retained,  it  should  be  retained  as  aTp  (eq.  (8b)).  In  the  problem  investi¬ 
gated,  p  turns  out  to  be  negligibly  small,  so  that  a  physically  significant 
error  would  develop  only  if  aT  »  1.  For  glasses  in  the  temperature  range 
of  the  problem,  aT  is  not  large,  so  that  no  substantial  numerical  errors 


result  from  this  cause. 
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